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^•Abstract 

The origin of Neptune's large, circular but retrograde satellite Triton has remained largely unexplained. There is an apparent 
jy-^ consensus that its origin lies in it being captured, but until recently no successful capture mechanism has been found. Agnor & 

Hamilton (2006) demonstrated that the disruption of a trans-Neptunian binary object which had Triton as a member, and which 
J~i underwent a very close encounter with Neptune, was an effective mechanism to capture Triton while its former partner continued 

"on a hyperbolic orbit. The subsequent evolution of Triton's post-capture orbit to its current one could have proceeded through 
1^ gravitational tides (Correia, 2009), during which time Triton was most likely semi-molten (McKinnon, 1984). However, to date, 

ho study has been performed that considered both the capture and the subsequent tidal evolution. Thus it is attempted here with the 
Onuse of numerical simulations. The study by Agnor & Hamilton (2006) is repeated in the framework of the Nice model (Tsiganis 
Q et al., 2005) to determine the post-capture orbit of Triton. After capture Triton is then subjected to tidal evolution using the model 
^ of Mignard (1979, 1980). The perturbations from the Sun and the figure of Neptune are included. The perturbations from the Sun 

Acting on Triton just after its capture cause it to spend a long time in its high-eccentricity phase, usually of the order of 10 Myr, 

1 ^ I while the typical time to circularise to its current orbit is some 200 Myr, consistent with earlier studies. The current orbit of Triton 

is consistent with an origin through binary capture and tidal evolution, even though the model prefers Triton to be closer to Neptune 
than it is today. The probability of capturing Triton in this manner is approximately 0.7%. Since the capture of Triton was at most 
^ a 50% event - since only Neptune has one, but Uranus does not - we deduce that in the primordial trans-Neptunian disc there were 
some 100 binaries with at least one Triton-sized member Morbidelli et al. (2009) concludes there were some 1 000 Triton-sized 
bodies in the trans-Neptunian proto-planetary disc, so the primordial binary fraction with at least one Triton-sized member is 10%. 
T— I This value is consistent with theoretical predictions, but at the low end. If Triton was captured at the same time as Neptune's 
ly-^ irregular satellites, the far majority of these, including Nereid, would be lost. This suggests either that Triton was captured on an 
orbit with a small semi-major axis a < 50 Rn (a rare event), or that it was captured before the dynamical instability of the Nice 
^-H tnodel, or that some other mechanism was at play. The issue of keeping the irregular satellites remains unresolved. 
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?H 1. Introduction and background belief that Triton was captured by Neptune from a heliocen- 

"^J trie orbit (McKinnon, 1984; Benner & McKinnon, 1995), and 

Of all the large natural satellites of the planets, Neptune's Neptune-centric orbit subsequently decayed through either 

Triton is one of the most pecuUar. It orbits Neptune at a dis- tidal interaction (McCord, 1966; McKinnon, 1984; Goldreich 

tance of 354 759 km or 14.3 Neptune radii (Rn), with a period of gt al., 1989; Correia, 2009), Neptune's circumplanetary gas disc 

5.877 days (Jacobson et al., 1991), similar to some of Uranus' (McKinnon & Leith, 1995) or a debris disc that formed from 

regular satellites. Its orbit is circular (e ~ 10-^) (Jacobson et collisions among Neptune's putative regular satellites (Cuk & 

al., 1991), yet its incUnation with respect to Neptune's equa- Gladman, 2005). We shall briefly discuss each of these scenar- 

tor is 156.8°, thus retrograde and ruling out a formation from ^elow. For an excellent in-depth review article on Triton 

within the Neptune system. Its retrograde orbit has led to the ^nd its origin we refer to McKinnon et al. (1995). 
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1.1. Capture 

Littleton (1936) hypothesised that Triton and Pluto origi- 
nated as adjacent prograde satellites of Neptune and that ejec- 
tion of the latter left the former on a retrograde orbit. Yet McK- 
innon (1984) demonstrated that this scenario is impossible: the 
amount of mass and angular momentum in the system is in- 
sufficient to make Triton retrograde. McKinnon (1984) argued 
instead that Triton and Pluto originated from a reservoir of bod- 
ies in the outer solar system. Triton was captured by Neptune 
during a close approach and its orbit was circularised through 
tidal evolution to its current one. McKinnon (1984) concluded 
that the tidal heating that followed as Triton circularised should 
have melted it. However, McKinnon's (1984) capture scenario 
presented a problem. In order to capture Triton it needs to lose 
enough energy to be bound from hyperbolic orbit and not be 
able to subsequently escape through the Hill sphere. This is 
difficult to do in one orbit using tidal interaction alone. There- 
fore Benner & McKinnon (1995) studied temporary capture 
from heliocentric orbit and the subsequent evolution of a set of 
these temporary-captured orbits in the circular restricted three- 
body problem consisting of the Sun, Neptune and a massless 
Triton. Solar perturbations acting on Triton cause its angular 
momentum to oscillate with a period half of Neptune's orbital 
period, with secular perturbations acting on longer time scales. 
In extreme cases these perturbations align and decrease Triton's 
pericentre distance, q, to within a few Neptune radii. Here the 
collision with an existing satellite or aerodynamic drag from a 
putative circumplanetary nebula could have dissipated enough 
energy to make capture permanent. However, Benner & McK- 
innon (1995) favour prompt capture at low q by collision or gas 
drag because a long-lasting temporary capture would have re- 
sulted in large changes in q and thus decreasing the probability 
of a collision or the eff'ect of gas drag. A prompt capture ensures 
there are many more close flybys which speed up the reduction 
of the orbit. In any case, permanent capture was difficult to 
achieve. 

A diff'erent scenario for the capture of Triton that did not 
require complex dynamics was proposed by Goldreich et al. 
(1989), who favoured the idea of a collision having occurred 
between Triton and a hypothetical regular satelUte of Neptune. 
Goldreich et al. (1989) argued that a collision with a regular 
satellite containing a few percent of Triton's mass would have 
shattered the satellite and left Triton bound to Neptune. How- 
ever, this scenario could be problematic. The collision prob- 
ability of Triton with a regular satellite is 10"^ per pericentre 
passage. In order to make this scenario work with a reasonable 
probability, Goldreich et al. (1989) argue there had to be 10 000 
Triton-mass objects encountering Neptune within 10 Rn, which 
is too large by several orders of magnitude compared to recent 
estimates derived from the size-distributions of the Kuiper Belt 
and Jupiter's Trojans (Morbidelli et al., 2009). However, it may 
not be necessary that each of these pericentre passages within 
10 Rn has to come from a unique object. Instead, there just 
have to be of the order of 10^ such passages of Triton-like ob- 
jects. The question then becomes whether this many passages 
is feasible. We shall return to this problem in the next section 



when we discuss our numerical simulations. Further criticism 
of this scenario came from McKinnon et al. (1995), who argue 
that unless the satellite that Triton would collide with is tiny and 
contains less than 2% of the mass of Triton, Triton would have 
shattered too. Since a collision with the larger satellites is much 
more likely, and assuming that Neptune's putative regular satel- 
lites were similar in size and mass to those of Uranus, neither 
the satellite nor Triton would have survived the collision. 

Given the difficulties of permanently capturing Triton in the 
above scenarios, Agnor & Hamilton (2006) suggested a new 
idea. They developed analytical arguments and used numeri- 
cal simulations to show that Triton could be captured through 
the dissociation of a binary planetesimal when it passed close 
to Neptune. These three-body encounters will disrupt the bi- 
nary if its centre of mass passes close enough to Neptune that 
the binary's orbital separation is approximately equal to its Hill 
sphere. Agnor & Hamilton (2006) tested this scenario with 
binaries consisting of objects with masses 1 nij and 0.1 nij, 
where nij is the mass of Triton. Agnor & Hamilton (2006) 
showed that capture is plausible for a large variety of initial 
conditions of the binary, such as the velocity at the time of en- 
counter with Neptune and the separation of the binary. It turned 
out that the probability to capture Triton decreased rapidly once 
the encounter velocity exceeded about 10% of Neptune's orbital 
velocity, while the lighter body was captured for encounter ve- 
locities of up to 70% of Neptune's orbital velocity. This result 
is not surprising because the change in velocity at disruption 
experienced by the heaviest body, Avi, is Avi oc m2, while for 
the lighter body Av2 <x nii (Agnor & Hamilton, 2006), and thus 
the fighter body experiences the larger kick and is more easily 
captured. It then automatically foUows that for fighter binaries, 
while more easily disrupted because of their lower binding en- 
ergy, it is also more difficult to capture one of its members: the 
low binding energy translates into a very small Av and one of 
them is bound to Neptune on a virtually parabolic orbit, where 
solar perturbations are stronger Therefore, capture prec eded by 
binary disruption is favoured for heavy binaries or those with a 
large mass ratio. Agnor & Hamilton (2006) found that the semi- 
major axis of Triton upon capture would always exceed 300 Rn, 
whUe the lighter body could be captured much closer to Nep- 
tune. Unfortunately, Agnor & Hamilton (2006) do not provide 
any statistics for their mechanism so that the probability of this 
event having occurred could not be tested nor compared to other 
results. 

The lack of statistical information about the success rate of 
capturing Triton through the disruption of a planetesimal binary 
led Vokrouhlicky et al. (2008) to investigate whether or not Tri- 
ton and the other irregular satellites of the giant planets could 
be captured via binary dissociation. This project was aimed at 
providing a comparison with the planet-planet encounter model 
of Nesvorny et al. (2007). The initial conditions of Vokrouh- 
licky et al. (2008) come from the Nice model (Tsiganis et al., 
2005), and they perform a series of simulations of the plane- 
tary instability, recording all tfie close encounters between tfie 
planets and planetesimals. They subsequently send a very large 
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number of binaries past each planet, with the distribution of 
the binary's hyperboHc planetocentric orbits taken from the en- 
counter parameters recorded earher, and their size distribution 
taken from that of Kuiper Beh objects (e.g. Bernstein et al., 
2004). They had difficulty capturing many small satellites from 
the dissociation of binaries, for reasons that we just explained 
above. When considering the capture of Triton they conclude 
that near-equal mass binaries produce the most captures, with a 
mass ratio of 1:2 to 1:3 preferred. While larger binary primary 
masses yield more captures, there are fewer of these around 
so that the overall capture probability decreases. Together with 
the fact that capturing members from a light binary is more diffi- 
cult, this could explain why we see no sub-Triton mass irregular 
satellites orbiting Neptune and Uranus. In conclusion, they find 
an overall capture probability of less than 2%, with no restric- 
tions on the semi-major axis, eccentricity or inclination. They 
add that Triton's capture was most likely to have occurred 5- 
10 Myr after Neptune's formation when the planetesimal disc 
was kept dynamically cold by the surrounding gas of the so- 
lar nebula; if it was captured during the planetary instability of 
the Nice model, it is also preferred that it is captured early on 
before the disc is too dynamically hot and thus the encounter 
velocity with Neptune is too large. 

Thus, it appears that the favoured mechanism for the capture 
of Triton is through the dissociation of a binary that had a very 
close encounter with Neptune. After capture, the orbit of Triton 
needs to shrink to its present size either through tides or other 
means. 

1.2. Post-capture evolution 

After Triton was captured and gravitationally bound to Nep- 
tune, several mechanisms have been invoked to evolve it to its 
current orbit. These are tidal interaction with Neptune (Mc- 
Cord, 1966; McKinnon, 1984; Goldreich et al., 1989; Correia, 
2009), interactions with a circumplanetary gas disc (McKinnon 
& Leith, 1995) or a debris disc (Cuk & Gladman, 2005). 

The first attempt at calculating Triton's orbital history with 
tidal interaction was performed by McCord (1966), who used 
the tidal model of MacDonald (1964) and expanded the equa- 
tions expressing Triton's change in orbital elements up to sixth 
order in eccentricity. By using reasonable estimates for Tri- 
ton's Love number (^2), quality factor {Q) and mass, McCord 
(1966) concludes that Triton could have reached its current or- 
bit from a highly extended, nearly parabolic orbit within the age 
of the solar system. The typical time to become circular is some 
100 Myr for Qj - 100, a typical value for rocky bodies. 

The results of McCord (1966) were verified by McKinnon 
(1984) and Goldreich et al. (1989). McKinnon (1984) argues 
that the tidal dissipation in Triton would have been enough 
to melt it, increasing the dissipation and thus shortening the 
time to reach its current orbit. Goldreich et al. (1989), on the 
other hand, use a simplified tidal model valid only for eccen- 
tricity e «c 1, but because of angular momentum conservation 
they argue that for e ~ I the tidal evolution only changes the 



semi-major axis a at constant pericentre, q. Assuming a value 
ki/Q = 10"^, typical for rocky bodies, they find that Triton 
could have reached its current orbit from a post-capture orbit 
with a - I 000 Rn within 400 Myr, but argue that a semi- 
molten Triton would have increased the rate of orbital decay. 
However, they added that the perturbations from the Sun cause 
repeated changes in the angular momentum with a period equal 
to half of Neptune's orbital period. These perturbations in an- 
gular momentum cause q to oscillate by as much as 10 Rn when 
a > 600 Rn- These oscillations greatly reduce the effect of tides 
and increase the time to circularise the orbit, but the authors do 
not specify by how much. 

Correia (2009) uses the tidal models of Mignard (1979, 1980) 
and Hut (1981), which do not require series expansions in the 
eccentricity, and he includes keeping track of Triton's obliq- 
uity and spin rate. By assuming k2/QT = lO^^* and an initial 
spin period of Triton of 24 h, Correia (2009) shows that Triton 
reaches its current orbit and spin rate in less than 1 Gyr from an 
orbit with a ~ 2000 Rn and q ~ 7 R^. In addition, he argues 
that tides are sufficient to account for Triton's low eccentricity 
(e ~ 10"^) and obliquity (st = 0°.46). Thus, from the above 
studies and with reasonable tidal parameters, it appears as if 
tides raised on Triton by Neptune are capable of reducing it to 
its current orbit from a post-capture one well within the age of 
the solar system, provided that Triton's pericentre upon capture 
is very close to Neptune (q ~ 7 Rn)- 

A different approach to reducing Triton's orbit from an ex- 
tended post-capture one was performed by McKinnon & Leith 
(1995), who studied the influence of a circum-Neptunian gas 
disc on the orbit of a post-captured Triton. They mimic a pu- 
tative nebula around Neptune after a minimum-mass Uranus 
nebula. They report strong evolution of the eccentricity and 
semi-major axis but almost no change in the inclination. Their 
results are insensitive to the radial surface distribution. Thus 
Triton could have evolved to its current orbit through gas drag 
and subsequent tidal evolution because, after the gas has disap- 
peared, Triton's eccentricity is about 0.2. Gas drag is able to 
reduce Triton's orbital angular momentum to its current level 
in about 1 000 years in the absence of solar perturbations. The 
latter can increase the gas drag time scale to lO'^-lO^ years, 
so that Triton could have outlived a hot, turbulent nebula last- 
ing some 1 000 years, but not a cool, long-lived low-mass one 
(10'' yr). In addition, McKinnon & Leith (1995) argue that Tri- 
ton could have cleared an annulus in the gas that could have 
halted its orbital decay. Thus, even though gas drag seems a 
favourable mechanism to circularise Triton from a post-capture 
orbit, it requires sensitive timing and the role of solar perturba- 
tions might reduce its impact. Recently, Ayliffe & Bate (2009) 
have performed the most sophisticated simulations of circum- 
planetary gas discs and they do not encounter circumplanetary 
discs around protoplanets with masses similar to Uranus and 
Neptune. This may not be a problem for Neptune because of 
the existence of Triton. However, Uranus' satellite system ex- 
hibits properties that are suggestive of a disc origin. The latter 
was made popular by Canup & Ward (2006) who discovered 
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that the typical mass ratio of the regular satellites relative to 
the giant planets should be ~ 10 The Uranian satellites are 
all approximately an order of magnitude less massive than that 
(Jacobson et al., 1992), suggesting a possible different origin 
than the one advocated by Canup & Ward (2006). One alterna- 
tive is that the satellites formed through the viscous spreading 
of an impact-generated debris disc (Ward & Canup, 2003). In 
another alternative Ayliffe & Bate (2009) suggest that a proto- 
planet's rotating envelope must cool following the dispersal of 
the encompassing circumstellar disc. As it does so it may well 
flatten into a disc suitable for satellite growth, suggesting a rela- 
tively late formation of the satellites. In any case, further study 
is needed to understand the formation of the Uranian satellite 
system and the idea that Triton was captured by gas drag. 

Both the tidal evolution model and the gas drag model have 
the disadvantage that they are most eff'ective very close to Nep- 
tune. However, it is likely that after capture Triton's semi -major 
axis was very large, a > 100 Rn- Beyond this distance, so- 
lar perturbations become important. Specifically, the Sun in- 
duces perturbations in Triton that are the same as the Kozai ef- 
fect experienced by high-inclination asteroids when perturbed 
by Jupiter (Kozai, 1962). The Kozai mechanism conserves the 
z-component of the angular momentum, L - cos /, 

where / is the orbital inclination of Triton with respect to the or- 
bit of the perturbing body (the Sun, or Neptune's orbital plane). 
Hence forprograde orbits inclination and eccentricity oscillate 
out of phase, while for retrograde orbits these oscillations are 
in phase. Depending on the initial conditions of the orbit, the 
Kozai effect can greatly increase the pericentre distance of Tri- 
ton after capture. In turn, this could greatly lengthen the time it 
takes for tides or gas drag to circularise the orbit. These prob- 
lems were pointed out by Cuk & Gladman (2005). They argued 
that the effect of the Kozai oscillations of Triton's pericentre in- 
creases the time scale to reach a circular orbit beyond the age 
of the solar system and therefore tides alone cannot circularise 
Triton. They reached their conclusion by averaging a simplified 
version of the tidal equations over one Kozai cycle. In order to 
find a different and faster mechanism to circularise Triton they 
investigate the role of an eccentric, retrograde Triton on Nep- 
tune's putative regular satellites. Cuk & Gladman (2005) argue 
that an eccentric Triton forces a high eccentricity on Neptune's 
putative regular satellites, which begin to cross each other and 
collide within a few thousand years. The collision would shat- 
ter both of these regular satellites and create a debris disc that 
Triton would pass through. By modelling the effect of the disc 
as a series of impulsive kicks occurring at pericentre, they con- 
clude that the time scale for the evolution of Triton is ~ 10^ yr, 
much shorter than for tides. The natural course of events is for 
Triton to sweep up all of the mass in the disc through collisions. 
The added benefit of this fast circularisation is that it could save 
Nereid from being ejected or colliding with Neptune. 

Thus, the debris disc of Cuk & Gladman (2005) seems an 
interesting alternative to the tidal or gas drag models. Unfor- 
tunately they do not specify in detail how they modelled the 



action of the debris disc on Triton so that its results cannot be 
verified. In addition, their arguments might not hold because 
of the following. When considering Uranus' satellite system, 
perturbations from Triton would increase their eccentricities to 
orbit-crossing values and their inclinations to a few degrees. 
Using the formulation of Opik (1976), the collision probabil- 
ity between any of the regular satellites is also 10"^ per orbit. 
Since the orbital period of the regular satellites is much shorter 
than that of Triton just after its capture, these satellites should 
suffer several collisions before any of them hits Triton. Most 
of these mutual collisions occur at impact velocities approxi- 
mately 1 to 2 times the escape velocity of the satellites. Agnor 
& Asphaug (2004) show that collisions with such a low mu- 
tual velocity tend to be either merging or hit-and-run collisions. 
For most of these hit-and-run collisions, the ejecta would be 
approximately 10% of the mass of the smallest satellite that 
is involved. Perturbations from the other satellites and Triton 
will quickly increase the eccentricities of the ejecta and most 
of these are rapidly swept up. It is likely that eventually a col- 
lision with Triton will occur that will shatter both the satellite 
and Triton. To summarise, it is not immediately clear whether 
these mutual collisions among the satellites will create the de- 
bris disc that Cuk & Gladman (2005) suggested, or whether the 
satellites will remain largely intact and the system is destroyed 
by a collision with Triton. 

In addition to a possible collision, there is another outcome: 
ejection of Triton. The forcing of the eccentricities and incli- 
nations of the regular satellites by Triton depends on both the 
eccentricity of Triton and the semi-major axis ratio between 
Triton and the satellites. For large semi-major axis ratios the 
eccentricities of the satellites could remain small enough to 
avoid crossing. Every time Triton passes through the satellite 
system it receives a root-mean-square energy kick of the order 
of A(l/a) ~ Gmr/a, (Fernandez, 1981; Duncan et al., 1987), 
where nir is the mass of Triton and a, is the semi-major axis 
of a satellite. The satellites receive an energy kick of the same 
magnitude. These kicks cause Triton to random walk in energy. 
Typically a, ~ 10 Rn and so A(l/fl) ~ 2 x 10"^ 7?^'. For a 
captured orbit similar to that of Nereid I /a ~ 4 x 10"^ R^^ and 
the number of orbits to random walk to ejection is then around 
40 000, ignoring a possible rare Levy flight, and thus the time 
to eject Triton is some 10^ years. These are crude estimates 
but what is important to notice is that the ejection time and the 
collision time are very similar. Since most collisions leave the 
majority of the mass in the satellites instead of debris, ejection 
becomes a feasible outcome. Clearly further study is needed 
to determine the most likely scenario when considering a pu- 
tative regular satellite system of Neptune being perturbed by a 
just-captured Triton. 

1.3. Current approach 

In this paper we investigate the capture and subsequent evo- 
lution of Triton in more detail, building on some of the earlier 
works. First, we study the capture of Triton using the binary 
capture scenario of Agnor & Hamilton (2006), but in the frame- 
work of the Nice model, as in Vokrouhlicky et al. (2009). This 
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model has booked a large number of successes, such as explain- 
ing the delay that caused the Late Heavy Bombardment of the 
terrestrial planets (Gomes et al., 2005), the origin of Jupiter's 
Trojans (Morbidelli et al., 2005), the structure of the Kuiperbelt 
(Levison et al., 2008b), the dichotomy of the asteroid belt (Lev- 
ison et al., 2009) and the irregular satellites of the giant planets 
(Nesvomy et al., 2007). Like Vokrouhlicky et al. (2009) we 
re-enact a series of encounters with binary planetesimals. This 
will give us the distribution of orbits of Triton just after capture. 
Second, we use the tidal equations of Correia (2009), based on 
the model of Mignard (1979, 1980), to determine which of the 
captured orbits are able to circularise within the age of the so- 
lar system, and which initial conditions will place Triton ap- 
proximately 14 Rn away from Neptune on a circular orbit. We 
improve upon Correia's (2009) tidal model, as well as those of 
McCord (1966) and Goldreich et al. (1989), by adding the per- 
turbations of the Kozai mechanism caused by the Sun, as sug- 
gested by Cuk & Gladman (2005), as well as the perturbations 
from Neptune's figure. The latter should suppress the Kozai ef- 
fect once the semi-major axis reaches ~ 73 Rn (Kinoshita 
& Nakai, 1991). The final outcome of the capture and tidal 
simulations should give us a probability to obtain a circular, 
retrograde Triton. This probability is then translated into the 
expected number of Triton-like objects in the proto-planetary 
disc for consistency. This paper is divided as follows: in sec- 
tion 2 we present our models and employed methods. Section 3 
contains the results from the numerical experiments. In sec- 
tion 4 we discuss how Triton's capture affects Neptune's irreg- 
ular satellites, concentrating on Nereid. Conclusions and sum- 
mary follow in the last section. 

2. Model and methods 

The model that we have employed consists of three parts. 
First, we record the encounters of planetesimals with Neptune 
during a Nice model simulation. Second, the deepest encoun- 
ters are re-enacted with binaries with varying mass ratio to de- 
termine the capture probability and the resulting orbital distri- 
bution of captured objects. Third, the orbits of the captured 
orbits are evolved under the action of tides to determine which 
ones end up being similar to Triton. Below we explain each 
stage in detail. 

2.1. Planet migration 

We consider a system consisting of the Sun, the four gi- 
ant planets and a planetesimal disc. The initial solar system 
was more compact than today, and the planets are thought to 
have formed between 5 and 15 AU on quasi-circular, coplanar 
orbits (Tsiganis et al., 2005). We took the initial conditions 
of Gomes et al. (2005) which ensured that Jupiter and Saturn 
crossed their 2:1 orbital resonance some 600 Myr after their 
formation, which subsequently triggered the planetary instabil- 
ity (Gomes et al., 2005). We assumed that the planetesimal 
disc was situated just beyond the orbits of the planets, ending at 
30 AU, with radial mass distribution varying as r"'. The plan- 
etesimal disc was composed of 10000 equal-mass bodies with 



35 




I ' ' ' 1 

0.7 0.75 0.8 0.85 0.9 

T [Gyr] 

Figure 1: Semi-major axis, perilielion and aphelion distance of tlie four giant 
planets around the instability. The letters designate which lines correspond to 
which planet. Neptune's path is traced with thick lines. 

a total mass equal to 35Me (where is the Earth's mass). 
We do not use the initial conditions that arise from the previ- 
ous phase of migration induced by the gas disc (Morbidelli et 
al., 2007) because we have a larger database of simulations of 
the former, and because the interactions between planets and 
planetesimals are similar in both discs after the late instabil- 
ity. We simulated the dynamical evolution using the numeri- 
cal integrator MERCURY (Chambers, 1999), where we use a 
computational 'trick' to decrease the amount of CPU time dur- 
ing the planetary migration simulations (Gomes et al. 2004). 
We defined an encounter to occur when the distance between 
a planetesimal and a planet is less than d - /Rh, where Rh is 
the planetary Hill radius and / is a factor larger than unity. We 
used d = I AU. Each encounter within this distance was regis- 
tered in detail, keeping track of the position and velocity of the 
planetesimal and the planet in the heliocentric reference frame, 
and the time of closest approach. Once a planetesimal entered 
the planet's Hill sphere, the system switched to planetocentric 
coordinates, where it was observed that, as expected, most of 
the encounters are hyperbolic. However, some of these encoun- 
ters are elliptical and therefore the planetesimal remains around 
the planet for a relatively long time during a temporary capture. 
The total integration time was 4.5 Gyr with a time step equal to 
0.4 yr. 

We performed several Nice model simulations and chose 
the one where the final orbits of the giant planets are closest to 
their current ones. Right from the beginning, some planetesi- 
mals from the disc encounter the outermost planet and conse- 
quently the planets slowly migrate, increasing their orbital sep- 
aration. This slow leakage from the disc continues until after 
788 Myr when Jupiter and Saturn cross the 2:1 resonance. This 
resonance crossing causes the system to destabilise and one of 
the ice giants (Neptune) is scattered into the disc, destabilis- 
ing it, and scattering planetesimals all over the solar system. 
In our preferred simulation, Saturn and the ice giants undergo 
many mutual encounters and the ice giants exchange orbits. The 
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Planet Number Total enc Elliptic Hyperbolic 



Jupiter 
Saturn 
Uranus 
Neptune 



5 866 

8 249 

9 891 
9 934 



216921 
288213 
1375 907 
1 748 867 



8 844 
4 874 
10997 
65 232 



208 077 
283 339 
1364910 
1683 635 



Table 1: First column: the total number of unique planetesimals that had close 
encounters with this planet. Second column: the total number of close encoun- 
ters between planetesimals and this planet. Third column: the total number of 
elliptical encounters between a planetesimal and this planet. Fourth column: the 
total number of hyperbolic encounters between a planetesimal and this planet. 



number of planetesimals having close encounters with the plan- 
ets decreased quickly after the instability, and their population 
decays approximately exponentially. Some 50 Myr after the 
instability the planets are in their actual configuration. Fig. [1] 
plots the evolution of the four giant planets around the time of 
the instability. The lines represent their semi-major axis, peri- 
helion and aphelion distance respectively. The letters designate 
which set of lines corresponds to which planet: Jupiter stays 
around 5 AU, Saturn around 9 AU, Uranus goes from 1 8 AU to 
20 AU and Neptune (outlined in thick lines) jumps from 12 AU 
to 30 AU. 

The encounter data from this preferred simulation is listed 
in Table [T] The first column lists the planet with which the 
encounters occurred. The second column lists the number of 
unique test particles, out of 10000, that suffered an encounter 
with this planet during the lifetime. The third column lists 
the total number of planetesimals that encountered the planet 
within 1 AU. The fourth column lists the number of encoun- 
ters that were elliptical while the last column lists the number 
of hyperbolic encounters. There are several things that should 
be pointed out. The first is that Uranus and Neptune undergo 
many more encounters than Jupiter and Saturn, because the 
former are scattered into the pristine planetesimal disc while 
Jupiter and Saturn are left behind. The number of unique parti- 
cles each planet encounters is similar to the ratios described in 
Fernandez (1997) based on planetesimals that were scattered 
by Neptune from the Kuiper belt. Specifically, Jupiter typi- 
cally ends up controlling 58% of the planetesimals, Saturn ap- 
proximately 22% (but it encounters 85% of them) and Uranus 
and Neptune approximately 10% each (so that each will en- 
counter about 90% of the planetesimals). A second feature is 
that Neptune encounters substantially more planetesimals than 
Uranus, partially because some particles that are under its con- 
trol never encounter Uranus and instead suffer many encounters 
with Neptune on their way to ejection, and partially because 
the proximity of Saturn to Uranus decreases the dynamical in- 
fluence of the latter Vokrouhlicky et al. (2008) also reported 
a substantially lower number of planetesimal encounters with 
Uranus than with Neptune. Last, the significantly larger number 
of elliptical encounters with Neptune vs Uranus is most likely 
the result of the higher encounter velocities of the planetesimals 
with Uranus than with Neptune. 
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Figure 2: Top panel: velocity distribution of particles that encounter Neptune 
after the instabiUty. The raw data is shown by the bullets while the solid line 
shows the best- fit Maxwellian with i',,, = 1.31 km s"'. Bottom panel: cumula- 
tive distribution in pericentre distance with respect to Neptune. 



For our purpose we are only interested in planetesimal en- 
counters with Neptune after the instability. For this study, we 
analyse the distribution of the planetocentric orbits of the plan- 
etesimals that had close encounters with Neptune only. Specif- 
ically, we are interested in the velocity distribution of the plan- 
etesimals as they encounter Neptune and the distribution of their 
closest approach distance to Neptune. It turns out that the veloc- 
ity distribution of the encounters is roughly Maxwellian, whose 
functional form is given by 



2 v2 
3-expl 



(4)- 



(1) 



with Vm the parameter velocity, which was found to be 1 . 3 1 km s " ' . 
In the top panel of Fig. |2] the bullets depict the velocity distri- 
bution of the planetesimals as they encounter Neptune. The 
best-fit Maxwellian is plotted as a solid line. The bottom panel 
of Fig. |2]plots the cumulative distribution of the peri-Neptune 
distance, q. For small distances the cumulative distribution is 
linear in q while beyond 200 Rn the fit scales as q^, suggest- 
ing that gravitational focusing is only important for very close 
encounters. Both plots are in good agreement with those pre- 
sented in Vokrouhlicky et al. (2009). 

The cumulative q distribution in the bottom panel of Fig. |2] 
can be used to constrain the collisional capture scenario of Gol- 
dreich et al. (1989), which was mentioned in the introduction. 
Given that the probability of Triton colliding with a fictitious 
Neptunian regular satellite is approximately 10"^ for each peri- 
Neptune passage, one needs of the order of 10^ passages of 
Triton-like objects within approximately 20 Rn in order to cap- 
ture it through a collision event. From the bottom panel of Fig.|2] 
the probability per encounter of a planetesimal coming within 
20 Rn is of the order of 2 x 10"^. Combined with the intrin- 
sic collision probability with these fictitious satellites Neptune 
needs to undergo of the order of 10^ encounters with Triton- 
like objects. Using the data from Table [1] as guidance, on av- 
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erage each planetesimal undegoes about 100 encounters with 
Neptune, so that for each Triton-Uke object the probabihty of 
colHsional capture is approximately 10"^. This low probabil- 
ity makes the collisional capture scenario seem a very unlikely 
mechanism to account for the existence of Triton. 

Now that we have a data base of encounters with Neptune, 
the next step is to re-enact these encounters using binaries and 
determine which binary systems are disrupted and also leave 
Triton bound to Neptune. 

2.2. Binary encounters 

After recording the number of encounters with Neptune from 
the migration simulation, we re-enacted these using binary sys- 
tems. The centre of mass of the binary coincides with the tra- 
jectory of a planetesimal that had a close encounter with Nep- 
tune. From Agnor & Hamilton (2006) we know that a binary 
is ionised once the separation between the members becomes 
equal to the Hill sphere of the binary. This tidal disruption dis- 
tance is given by 



rtd 
Rn 



3pwW nil 
pi j\mi +m2 



1/3 



Ri' 



(2) 



where Rf^ is Neptune's radius, Ri is the radius of the primary 
binary component (usually Triton's radius), is the mean den- 
sity of Neptune, pi is the mean density of the primary compo- 
nent of the binary (usually Triton's density) and nii and nii are 
the masses of the primary and secondary of the binary. The 
approximation sign in equation ^ above is valid when consid- 
ering densities appropriate for Neptune and Triton. We created 
four different groups of binaries with different mass ratios, with 
each binary containing one member with a mass equal to Tri- 
ton's mass. The other member was either 0.1 nir, 0.3 nij, 1 nix 
or 3 niT- All binaries had an initially circular orbit and semi- 
major axis flg = 1.5 Rn = 37 146 km, which is the approxi- 
mate maximum value for the currently-known trans-Neptunian 
binary population (Noll et al., 2008). These binaries should 
disrupt once they come closer to Neptune than q ~ 27 R^- The 
other orbital elements were chosen at random. For each binary 
group we simulated 1 000 different orbits. To be on the safe 
side, we re-enacted only encounters with q < 100 R^. 

After we created the binary system, we simulated their ap- 
proach with Neptune using the MERCURY integrator, and reg- 
istered the changes in semi-major axis, eccentricity and inclina- 
tion of the binary. We consider that a member is ejected from 
the system whenever it has a > 3000 Rn- If the pericentre dis- 
tance is smaller than the planetary radius the body is removed 
from the simulation. 



2.3. Tidal evolution 

We decided to use the tidal equations presented in Correia 
(2009), which are based on the tidal models of Mignard (1979, 
1980) and Hut (1981). This model assumes that the time delay 
of the tidal response of a body. At, is a constant. This model 
can be made linear in At and the tidal equations can be written 
in closed form for all values of eccentricity e < I. However, 
the constant time delay model does not appear to agree with 
geophysical data (Efroimsky & Lainey, 2007). The rate of en- 
ergy dissipation within a body is characterised by the tidal pa- 
rameter Q, which is a measure of how many tidal oscillations 
are needed to damp the energy by order of itself. The constant 
time delay model of Mignard (1979, 1980) assumes that Q is 
inversely proportional to the frequency with which the body is 
distorted, x- For eccentric orbits x is the angular frequency of 
this body at pericentre. Laboratory experiments indicate that 
Q cc x", where a ~ \ for the range of periods that we are in- 
terested in (Karato, 2008). Nevertheless, we are mostly dealing 
with orbits with e ~ 1 where series expansions are invalid, and 
thus we decided to adopt this model because it best suits our 
needs. We decided to use a simplified model in which we inte- 
grate the equations governing the changes in semi-major axis, 
eccentricity, spin rate and obliquity of Triton, while ignoring 
any changes in the rotation rate of Neptune and its obliquity. 
This is justified because Correia (2009) found that the change 
in both of these quantities is negligible. We found that tides 
raised by Triton on Neptune are a couple of orders of magni- 
tude weaker than tides raised by Neptune on Triton, but we in- 
cluded them for the sake of completeness. The tidal equations 
are given by (Correia, 2009) 
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where sr is Triton's obliquity with respect to its own orbit, ir 
Triton's inclination with respect to Neptune's equator, a>r is Tri- 
ton's spin rate, cl>n is Neptune's spin rate, nr is Triton's mean 
motion, Cr is Triton's moment of inertia along its spin axis, and 



This part of the model allowed us to create a distribution 
of orbits that Triton could have had just after capture. The 
next step is to evolve these orbits using tides raised by Nep- 
tune on Triton to determine what fraction of post-capture orbits 
will yield the current orbit of Triton within the age of the solar 
system. 
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Figure 3: Tidal evolution of semi-major axis and pericentre distance of a Triton- 
sized object captured by Neptune. Initial conditions are from Correia (2009). 
Solid lines depict the evolution from Con'eia (2009) while the dashed lines 
show our simplified implementation. Solid lines data set courtesy of Alexandre 
Con'eia. 

= (l+3e2 + 3eV8)(l-eV^^ 
hie) = (l-t-15eV2 + 45e''/8 + 5eVl6)(l-e^^^ 
/3(e) = (l-i-31eV2 + 255eV8-i-185£'Vl6 + 25eV64) 

/4(e) = (1+3^2/2 + ^4/8X1-^2)-', 

/5(e) = (l + 15e2/4+15eV8 + 5eV64)(l-e2)-i3/2_ 

Here krr is Triton's Love number ^2 and b^tj is Triton's tidal 
response time, Rj is Triton's radius and mj is its mass. The 
quantities carrying subscripts N are for Neptune. Table |2]lists 
the values of the various quantities used above. Most are taken 
from Correia (2009), and references therein. The values of ujj 
and sj are starting values, which evolve to their current ones. 
An example of the tidal evolution in this simplified model is de- 
picted in Fig. [3] and shows the evolution of semi-major axis and 
pericentre distance of a Triton-sized object captured by Nep- 
tune. The initial conditions are taken from Correia (2009) and 
are qq - 2354 Rf^, qo - 7 Rn, k - 157° and w,o - 1/24 hr"'. 
The solid lines show the evolution depicted in Correia (2009) 
while the dashed lines show the evolution according to our sim- 
plified model. The slight differences are caused by Correia 
(2009) taking Cassini states and properly taking Triton's obliq- 
uity evolution into account. 



Quantity Value 



niT 


1.0837 X 10"** Mo 


rriN 


5.156 X 10"^ Mq 


kiT 


0.1 


kiN 


0.407 


Rj 


9.044 X 10"'' AU (1353 km) 


Rn 


1.6554 X 10""* AU (24764 km) 


Ct 


Q35mTR\ 


(jjj 


6884.65 rad yr"' {Pj = 8 h) 


a>M 


3418.82 radyr-' (Pn = 16.11 h) 


St 


170° 


En 


28.56° 




1.02 s(eA, =9000) 


Atr 


808 s or 8080 s 


a& 


30.1 AU 




0.008 


Ji 


3.343 X 10-3 



Table 2: Values of various quantities that enter the tidal equations, as well as 
those of the Kozai mechanism and figure of Neptune. 

to the Moon and Mars. For the Moon, the inferred value of ^2 
from laser ranging is ki - 0.02664 and Q ~ 30 (Williams et al., 
2005). This low value of Q is caused by tidal dissipation within 
the Moon's core, and results in ^2/G ~ 0.0011. For Mars, its 
value of k2 ranges from 0.11 to 0.16 (Marty et al., 2009) while 
its value of Q ranges between 80 (Lainey et al., 2007) and 92 
(Yoder et al., 2003), suggesting that for Mars k2/Q ~ 0.0015. 
Thus, it seems that for rocky bodies ko/Q ~ 0.0013 while for 
semi-molten bodies ^2/G ~ 0.015. We decided to adopt an 
intermediate value of ^2 = 0.1 for Triton (McKinnon et al., 
1995) and change the value of Qt by a factor ten between 
the solid and semi-molten state. For the solid state we adopt 
ki/Qr - lO^^*, implying Qt - 100, and for the semi-molten 
state kilQT - IQ'^, implying Qt - 10. This corresponds to a 
time delay At - 808 s in its current orbit for the solid state, and 
8080 s when semi-molten. 

As stated in the introduction, the tides are not acting alone. 
For large orbits, the Kozai mechanism (Kozai, 1962) is at work, 
caused by perturbations from the Sun, and we need to take its 
influence into account (Cuk & Gladman, 2005). The Kozai 
mechanism has two constants of motion (Kozai, 1962; Kinoshita 
& Nakai, 2007) 



We considered two values of Af 7". The first, AIt = (QrWr)"' 
808 s, where we used the current value of ojt and corresponds 
to imposing Qt - 100, a commonly adopted value for rocky 
bodies and close to the inferred value for Mars of Qm - 92 
(Yoder et al., 2003). However, Triton was most likely semi- 
molten during its tidal circularisation (McKinnon, 1984), which 
decreases the value of Q and thus increases that of Afj-. Using 
lo as an example, it has kolQi ~ 0.015 and thus Qi ~ 3 for ^2 = 
0.05 (Lainey et al., 2009) but when using the fluid Love num- 
ber, appropriate for semi-molten bodies, k2 - 1 .292 (Anderson 
et al., 2001) and we have Qi ~ 86. For solid bodies, we can turn 



Hk = y[(2 + 3e2)(3cos-/- 1) + 15e2sin2/cos2w], 
h, = (l-e^y^cosl, (5) 

where we omit the subscript T since we are only dealing with 
Triton. Here Hk is the averaged Kozai Hamiltonian (Kozai, 
1962), li~ is the z-component of the orbital angular momentum, 
/ is Triton's inclination with respect to Neptune's orbital plane, 
CO is Triton's argument of pericentre and 



16(1-4)3/2- 



(6) 



The Kozai mechanism induces coupled oscillations in eccen- 
tricity and inclination, and either a circulation or libration of co, 
depending on the value of Hk and h~. The equations of motion 
of these three variables, and Q, the longitude of Triton's ascend- 
ing node on Neptune's orbital plane (see below) are (Kinoshita 
& Nakai, 2007) 



to Neptune, we stopped integrating the Kozai and J2 effects be- 
cause the former is suppressed by the J2 precession and the 
latter only causes circulation of O and u), which do not con- 
tribute to the tidal evolution. Switching off the integration of 
these quantities significantly sped up the final part of the tidal 
evolution. 
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In addition to the Kozai mechanism, the figure of Neptune 
perturbs Triton's orbit once it is close enough to Neptune. While 
the perturbations of Neptune's figure do not alter the eccentric- 
ity and inclination, it does force a precession of the argument 
of pericentre. Once Triton is close enough to Neptune, this pre- 
cession will overtake the precession induced by the Kozai effect 
and the coupled oscillations in inclination and eccentricity will 
stop. Thus, in order to model the disappearance of the Kozai 
mechanism, we need to add the perturbations of Neptune's fig- 
ure to the tidal model as well. The averaged Hamiltonian is 
(e.g. Kinoshita & Nakai, 1991) 



Hn — — 



1 Gmn /RNfOcos^i- 1) 
A~) (l-e2)3/2 ' 
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where ; is Triton's incUnation with respect to Neptune's equator 
and J2 is Neptune's quadrupole moment. In order to add this to 
the tidal equations of motion and those of the Kozai mechanism, 
we need to transform Triton's inclination with respect to Nep- 
tune's equator to its inclination with respect to Neptune's orbital 
plane. This is done via cos / = cos cos / - sin e^r sin / cos Q, 
where is Neptune's obliquity. Unfortunately this adds an- 
other variable, O, to be integrated, and therefore the Kozai part 
of the regression of this angle needs to be taken into account as 
well (hence its inclusion above). We have 
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We integrated the system of equations consisting of (|3]l, (|7|i 
and ^ using a Bulirsch-Stoer integrator with variable time step 
(Bulirsch & Stoer, 1966). Once Triton was closer than 20 Rn 



When considering pure tidal motion, the orbital angular mo- 
mentum of Triton is approximately conserved. Thus, for ini- 
tially very eccentric orbits, the final, circular, orbit is located at 
a X 2qo, where qo is the original pericentre distance. Triton's 
current orbit suggests it was captured with go ~ 7 Rn- How- 
ever, the Kozai mechanism induces oscillations in the eccen- 
tricity and thus in q. From equations (O we have d oc q^^^/^ + 
(9(^-13/2-)^ so that a typical time scale on which the tides act, 
Ta = a/d cc q^^/^+0{q^^^^) is a very steep function of q and sug- 
gests the annulus in which the tides can circularise Triton within 
the age of the solar system is very narrow. Indeed, increasing 
Ta by an order of magnitude requires a relative increase in q 
of only 35%. Experimentation showed that the time to become 
circular reaches the age of the solar system for Qt - 10 when 
qo = 20 Rn, and thus orbits for which q never dips below 20 Rn 
can be ignored. Given the rapid increase in Ta with q, in the 
simplest and crudest sense one can envision the q dependence 
as a step function, where the tides are switched off if q is larger 
than some threshold value, and the tides are active when q is 
smaller. This behaviour suggests that the final orbit of Triton 
would have a ^ 2q^i„, where q^in is the minimum value of q 
that is obtained during the Kozai cycle. Indeed this appears to 
be a better approximation than the previous one a ^ 2qo. 

Now that we have all the ingredients in place, we report the 
results of our experiments below. In what follows, we removed 
any objects that achieved a final semi-major axis a < 5 Rn since 
it would then collide with Proteus (located at 4.8 Rn), or which 
encountered q < I Rn during their Kozai cycle, or whose q 
never reached below 25 Rn within 1 Myr (by comparison: the 
Kozai cycle time is typically 0.01 Myr). While we did inte- 
grate cases with Qt = 100, we do not present the results here 
since it was most likely that Triton was semi-molten (McKin- 
non, 1984). In what follows, the time scale for the tidal evo- 
lution of Triton should be considered as indicative rather than 
absolute. This is caused by the lack of knowledge of the tidal 
parameters and the model's simplification of the real tidal evo- 
lution. 

3. Results 

In this section the results from our numerical simulations 
are presented. First we present the distribution of captured or- 
bits. This is followed by a case study of the tidal evolution, after 
which we turn to the tidal evolution of all the captured orbits. 

3.1. Captured orbits 

Fig. |4] displays the pericentre distance (q) vs. semi-major 
axis (a) of objects with Triton's mass that were successfully 
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Figure 4: Pericentre distance vs. semi-major axis of captured satellites with 
Triton's mass after the disruption of a binary. These are orbits only upon capture 
without subsequent tidal evolution. 



Figure 5: Cumulative distributions of semi-major axis (top-left), pericentre 
(top-right), inchnation (bottom-left) and eccentricity (bottom-right) for Triton- 
mass satellites upon capture. 



captured following the disruption of a binary, as outlined in 
Section 2.2. Since the binaries were enacted from encounters 
with Neptune that had q < 100 Rn, it is no surprise that most 
captured objects have a rather small value of q. The truncation 
of the semi-major axis at approximately 3000 is caused by 
our condition that objects with a > 3000 Rn are considered un- 
bound since they will be outside Neptune's Hill sphere at apoc- 
entre. Fig.|5]shows the cumulative distributions of semi-major 
axis (top-left), pericentre (top-right), inclination (bottom-left) 
and eccentricity (bottom-right) for Triton-mass satellites upon 
capture. The steep increase in the distribution of q between 
8 Rn and 32 Rn is probably caused by our choice of binaries. 
Similarly, the sharp rise in the eccentricity distribution when 
e > 0.9 is might also be an artefact of our initial conditions. 
The inclination distribution shows there is a clear preference 
for retrograde orbits upon capture, probably because of the in- 
creased stability of retrograde orbits with respect to the size 
of Neptune's Hill sphere (Hamilton & Krivov, 1997). In sum- 
mary, most orbits are captured with a 6 (100,2000) Rn and 
q e (4, 32) Rn with a clear retrograde preference. In contrast, 
Agnor & Hamilton (2006) find that the median semi-major axis 
of Triton after capture is a > 1 000 Rn, while we find a value 
of fl ~ 500 Rn, suggesting that the encounter parameters when 
Neptune migrates are different from the static case. These orbits 
now need to be evolved using the tidal model presented above. 
We define a successful case when Triton achieves e ~ 10"^ in 
less than 4 Gyr In principle one could define a successful case 
as any capture for which e < 1 and q> \ Rn but then the current 
orbit is not reproduced and thus we decided to use the former 
criterion. 

3.2. Tidal evolution: Case study 

In order to determine whether or not we can obtain Triton's 
current orbit within the age of the solar system from the cap- 
tured orbits presented in the previous subsection, we need to run 
them through the tidal equations. In this subsection we show a 
representative case of the evolution of Triton upon capture. In 
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Figure 6: Example of tidal evolution of one object using two different models. 
The top panel shows the evolution in semi-major axis. The solid Hne uses 
Neptune's cuirent obliquity. The dashed line has Neptune's obliquity set to 
zero. The bottom panels show the evolution of the pericentre. Left panel has 
the current obliquity of Neptune while the right panel has Neptune's obliquity 
set to zero. 

the next subsection we present the results of the final orbits of 
all captured objects. We have decided to model the tidal evolu- 
tion using two different methods, in order to compare them. For 
the first method, we integrate the tidal equations from section 
2.2 using as the initial conditions the captured orbits presented 
above. We included Kozai and Neptune's figure and set Nep- 
tune's obliquity equal to its current value. With this method 
one can obtain Triton's final inclination with respect to Nep- 
tune's current equator For the second case, we use the same 
equations as above but set Neptune's obliquity equal to zero. 
This was done to determine how the Kozai perturbations from 
the Sun determine the final outcome. 

Fig lUpresents the results of the tidal evolution of one object 
using Qt - 10 i.e. the semi-molten case. The starting condi- 
tions are a = 648.2 Rn, q = 8.38 Rn, I = 52.1°, to = 207.1°. 
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The top panel shows the evolution of the semi-major axis vs. 
time. The solid line depicts the case with Neptune's current 
obliquity = 28.56° (case 1). The dashed line depicts the 
case where the obliquity of Neptune is zero (case 2). In both 
cases the pericentre distance oscillates with large amplitude and 
short period compared to the tidal evolution. The latter occurs 
quickest when is at a minimum, and when this happens the 
semi-major axis shrinks. With each decrease in semi-major axis 
at minimum q, the value of li~ decreases and the excursions in 
q (and e) are less pronounced. Eventually the semi-major axis 
shrinks to the point where the perturbations from Neptune's fig- 
ure compete with the Kozai cycles induced by the Sun. This oc- 
curs at approximately 70 Rn- Once the semi-major axis shrinks 
below this value the Kozai cycles cease and the rest of the evo- 
lution proceeds at constant angular momentum. However, this 
value of the angular momentum is not equal to the value upon 
capture. The final semi-major axis is approximately Uf ~ 2q^i„. 

Fig.|7]presents a zoom for case 1 around the time when the 
Kozai cycles stop. The dotted line plots the semi-major axis vs 
time while the solid line represents the pericentre distance, q. 
The oscillations in q decrease in ampHtude and the minima are 
consistently closer to Neptune. The decrease in the minimum 
distance of q is caused by the perturbations from Neptune's fig- 
ure not being aligned with those of the Sun. Far from Neptune 
the orbital angular momentum vector of Triton precesses per- 
pendicular to Neptune's orbital plane. As the semi-major axis 
of Triton shrinks, the perturbations from Neptune's figure be- 
come stronger compared to the Kozai cycles induced by the 
Sun. Neptune's equator is not aligned with its orbit and the 
perturbations from Neptune's figure force Triton's orbital an- 
gular momentum to precess perpendicular to Neptune's equator 
rather than its orbit. The transition happens around 70 Rn- As 
the angular momentum gradually starts to precess about Nep- 
tune's rotational pole, the value of h~, which is a constant for the 
Kozai motion, has to decrease, causing an increase in the maxi- 
mum eccentricity and thus a decrease in the minimum value of 
q. Eventually the minimum q has dropped close to 4 R^ and the 
tidal evolution is then so rapid that the semi-major axis shrinks 
considerably before Kozai cycles lift the pericentre again. How- 
ever, by this time, the amplitude of the Kozai cycles have de- 
creased considerably so that future cycles are quickly damped 
and the semi-major axis continues to decrease. By now the tidal 
evolution more or less conserves the current total angular mo- 
mentum since the Kozai cycles have stopped, and the orbit cir- 
cularises at fl ~ 2^,nin. 

The outcome for the two simulations are different and are 
presented in Table |3] The first column represents the elements 
fl, q and time to become circular (Tc). The second column lists 
the initial values. The third column marks the first case (current 
obliquity of Neptune). The fourth column marks the second 
case (no obliquity of Neptune). As one can see, the first case 
yields an orbit much closer to Neptune than the second case. 
The time to become circular varies by a factor of two. The dif- 
ference between the final semi-major axis and values of Tc are 
the result of the difference in Neptune's obliquity. For the first 
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Figure 7: Zoom of the tidal evolution of Fig.|6]around the time the Kozai mech- 
anism stops. The dashed line shows the semi-major axis vs time while the solid 
line shows the pericentre distance. The minimum q decreases, due to a mis- 
alignment between Neptune's rotational pole and orbital pole, until eventually 
it reaches below 4 R^- The tidal evolution is then so rapid that the semi-major 
axis shiinks sufficiently for the Kozai cycles to stop. The remaining evolution 
proceeds at constant angular momentum. 

Element Initial 1 2 



a 648.2 9.03 15.68 
q 8.38 9.03 15.68 
Tc 242.3 631.8 

Table 3: The initial and final semi-major axis (a), pericentre (q) and time to 
circularise {T^) for the sample orbit with two different tidal models. The column 

(1) depicts the fist case where Neptune has its current obliquity. The column 

(2) coiTesponds to the case where Neptune's obliquity is zero. 

case the orbital evolution proceeds very quickly once q~ \ R^, 
which does not occur in the second case, where q never drops 
below ~ 7 R^. 

Now that we have given an overview of the tidal evolution, 
we turn to what the final solutions are when the ensemble of 
captured orbits are run through the tidal model. 

3.3. Tidal evolution: Final orbits 

In this subsection we shall present the results of the possible 
final orbits of Triton after tidal evolution of the captured orbits. 
We shall focus on those cases that reach e - 10"^ within the 
age of the solar system. Other cases are discarded since they 
are incompatible with Triton's current orbit. 

Fig.[8]shows the cumulative distribution in semi-major axis 
of those objects that reach Triton's current eccentricity within 
the age of the solar system for the semi-molten case (Qt = 10). 
The solid line is for case 1, the dashed line represents case 2 (in- 
dicated on the plot). There are two important features to note 
here. First, case 1 yields many orbits with a small final semi- 
major axis because of the feature displayed in Fig.Qabove: the 
minimum value of q decreases as a decreases and the orbits 
circularise at approximately 2^niin- The second feature is that 
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Figure 8: Cumulative distribution of final semi-major axis for objects which 
reach e = 10"^ within the age of the solai' system. The solid line is for case 1, 
the dashed line represents case 2. 
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Figure 9: Top panel: Cumulative distribution of time to become circular. The 
solid fine is for case 1, the dashed line represents case 2. Bottom panel: time to 
become circular' vs. final semi-major axis. Bullets present case 1, open circles 
are case 2. The border between the filled and open spaces of the plot scales 
approximately as Of oc r^''^. 

Triton's current orbit at 14.3 Rn is always in the upper quartile 
of the distribution; in the worst case it is in the upper 5% of 
the distribution. This result would suggest that if tides were the 
dominant mechanism behind circularising Triton's orbit after 
capture, one would think that Triton should be closer to Nep- 
tune than it is today. 

Fig.|9]depicts the cumulative distribution of the time it takes 
for Triton to become circular (top panel) and the final semi- 
major axis vs. circularisation time in the bottom panel. Once 
again in the top panel the solid line presents case 1, the dashed 
line is case 2. As one can see, case 1 has a much shorter circu- 
larisation time than case 2. This is no surprise because the cir- 
cularisation time depends heavily on q as oc q^^^^ and since 
case 1 experiences the lowest values in q it should circularise 
the quickest. From the plot one can infer that the median time 



to become circular ranges from 30 Myr to 70 Myr depending on 
the tidal model that is used. However, care has to be taken when 
using this number because it does not specify the final value of 
semi-major axis that is obtained. For that we turn to the bottom 
panel, which depicts the final semi-major axis vs. the circulari- 
sation time. The bullets present case 1, the open circles are case 
2. The line between the filled and open parts of the panel scales 

2/15 

as fl/ oc Tj,' . For Triton at 14 R^, the time to become circular 
by tides ranges from approximately 200 Myr for most orbits, up 
to 4 Gyr for a few exceptional cases, depending on the configu- 
ration of the original orbit. Thus, Triton would have reached its 
current orbit approximately 200 Myr after its capture. This time 
scale is consistent with earlier results of McCord (1966), McK- 
innon (1984), Goldreich et al. (1989), McKinnon et al (1995) 
and Correia (2009), who all give a typical time scale of 100- 
500 Myr for Triton to become circular, even though not all of 
the above works use the same tidal model as ours, nor do they 
take the Kozai mechanism into account. However, our tidal cir- 
cularisation time scale is different from that predicted by Cuk & 
Gladman (2005), even if it agrees with other estimates in the lit- 
erature. We can think of two reasons for this discrepancy. The 
first is that Cuk & Gladman (2005) use a very simplified model 
for the tides acting on Triton, which does not have the same 
scaling with q as our formulation. Their tidal damping time 
scale proceeds as proportional to q^ while ours goes as q^^^^. 
Secondly, Cuk & Gladman (2005) average the tidal equation 
over one Kozai cycle. While from their description it is not 
clear how this is done, one would suspect that this results in us- 
ing an averaged value of q in the tidal equation rather than the 
minimum value. The latter controls the time scale for circulari- 
sation. 

After its capture Triton undergoes at least another 3.6 Gyr of 
tidal evolution due to tides raised on Neptune by Triton. Even 
though the influence of tides raised on Neptune by Triton are a 
couple of orders of magnitude weaker than tides raised on Tri- 
ton by Neptune, the subsequent evolution after Triton becomes 
circular cannot be ignored. For a circular orbit of Triton raising 
tides on Neptune we have 



2Ki 

m, 



Kn I cos iraJN _ A 



(10) 



Since cosijojM lyij ^ 1, we can ignore the second factor in 
the brackets, which makes the equation integrable with solution 
Qf - (aJ^^^H-Cf)^^'-' where C = 6k2NinTRl/^tNCL)N cos ij ylGjni^. 
Note that for retrograde orbits C < and the orbit shrinks. Over 
3.6 Gyr Triton's orbit should have shrunk by about 1 .3 Rn, from 
15.6 Rn to its current value, making Triton's final position com- 
patible with the tidal model at the 20% level or less. 

Finally, Fig. [TO] shows histograms of the distribution of the 
final inclination with respect to Neptune's equator of all circu- 
lar orbits that used case 1. There is a large surplus of retro- 
grade objects and the majority of these are situated between 
140° and 150°, slightly lower than Triton's current value of 
157°. The large number of objects in this bin can partially be 



12 



180 




Figure 10: Histograms of tlie frequency of tlie final inclination with respect to 
Neptime's equator for case 1 i.e. with Neptune's current obliquity. 

explained by the following. The inclination of Triton with re- 
spect to Neptune's equator (/) is related to its inclination with 
respect to Neptune's orbit (/), Neptune's obliquity (sn) and the 
longitude of the ascending node (Q.) by cos i = cos e^v cos / - 
sin Sat sin/ cos O. Since cos / is almost a constant for bodies 
close to Neptune, the value of / varies between ; + and 
i - sn with the circulation time of Q. Consequently, this same 
final value of / can arise from orbits with original inclinations 
I + sn and I - Bf^, where the final inclination depends on the 
original value of Q.. The maximum value of / is 180° and so 
180° - sw = 151.5°. Thus, the bin between 140° and 150° can 
sample orbits from the full range (/ + e^,/ - en), but orbits 
with higher / cannot. Alternatively, orbits with lower values 
of ; sample orbits with / < 120°, where the Kozai mechanism 
operates strongly and can drive orbits to collide with Neptune. 
These collisions are the reason for the paucity of orbits with fi- 
nal ; ~ 90°, even if these orbits are protected from Kozai mech- 
anism at small semi-major axis because of Neptune's Jo- Sim- 
ilarly, the maximum between 30° and 40° can be explained in 
a similar manner, though the number difference with the corre- 
sponding retrograde case is most likely a result of the retrograde 
cases being stable up to larger distances from Neptune. 

We close this subsection by presenting some statistics for 
the case where Neptune's obliquity is not zero since this is the 
most probable. Of all the Tritons that we analysed, 3% are cap- 
tured with an orbit having q < \ R^, so that these collide with 
Neptune promptly after their capture. Of the remaining pop- 
ulation there are three possible outcomes; circular with final 
semi-major axis a > 5 Rn, eccentric orbit with final semi-major 
axis fl > 5 Rn, and tidal evolution to inside of 5 Rn before be- 
coming circular The last category consists of approximately 
33% of the total population, so that some 67% survive on orbits 
with final semi-major axis a > 5 Rn- This surviving population 
has only 29% on orbits with e > 10"^ or 19% of the total pop- 
ulation. Combining these numbers it turns out that a circular 
Triton with final a > 5 Rn occurs approximately 50% of the 
time, with twice as many retrograde vs prograde cases. If Nep- 



tune's obliquity were zero, the above statistics are virtually the 
same. 

3.4. Probabilities 

Next we examine whether or not the existence of Triton 
is compatible with a capture during the planetary instability. 
Here we just concern ourselves with how many bodies of Tri- 
ton's size are inferred to have existed from Triton's current or- 
bit. However, before we continue we should pause for a mo- 
ment and reflect on what we are interested in. The question 
is whether or not the capture of a Triton-like body is plausi- 
ble, but not be so common that Uranus would have experienced 
the same event. We can make this argument even stronger: the 
Uranian satellite appears to be unperturbed and regular while 
Neptune's system appears to have been disrupted. For this to 
occur it is unimportant if the perturbing body was captured pro- 
grade or retrograde, and whether or not the orbit is circular at 
the current epoch. Only if the capture is a probable event is it 
worthy to ask the question about the most likely end states. 

The probability of Triton's capture through the dissociation 
of a binary, and its subsequent evolution, is split into the prod- 
uct of two probabilities: the probabiUty that over its whole life- 
time a binary suffers an encounter with Neptune close enough 
for it to be disrupted (Pi), and the probability that the binary 
member of Triton's mass is actually captured (Pa)- Their prod- 
uct is Triton's capture probability. We could restrict ourselves 
to only circular and retrograde outcomes, and thus we need to 
multiply the above total probability by the probability that after 
capture Triton reaches a circular orbit at the current epoch and 
is retrograde (P3). However being circular and retrograde at the 
current epoch is not necessary for disrupting Neptune's primor- 
dial satellite system. 

The probability Pi should reflect the encounter history of 
the binary until its distruption. Each encounter with Neptune 
softens the binary and, provided it is not disrupted, after many 
encounters it has no knowledge of its original binding energy 
(Parker & Kaverlaars, 2010). However, we cannot take this 
history into account in our current methods but given that the 
probability of passing Neptune close enough to be disrupted is 
approximately 10 and that each planetesimal undergoes ap- 
proximately 100 encounters, we believe that our method is ac- 
curate enough. It is not necessary to disrupt the binary on its 
first passage and thus we have used the total probability of dis- 
ruption averaged over many encounters rather than just for the 
first encounter 

The value of Pi depends on the distance from Neptune at 
which the binary gets disrupted, which in turn depends on the 
mass ratio and total mass of the binary. Similarly, P2 depends 
on the system configuration to determine how many Tritons are 
captured. In table|4]we list the values of Pi and P2 for the binary 
systems that we considered, in units of percent. The values of 
Pi are for encounters that occur during and after the planetary 
instability, since we have information about how much mass 
there was available at this time (Gomes et al., 2005). Since the 
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Masses [wj-] m [R^] Pi [%] Po [%] ^"1 x P2 [%] 



1.0, 0.1 


35 


3.9 


19.3 


0.75 


1.0, 0.3 


33 


3.7 


19.8 


0.73 


1.0, 1.0 


29 


3.1 


36.9 


1.1 


3.0, 1.0 


23 


2.1 


11.8 


0.25 



Table 4: Table with probabilities Pi and P2 and the total capture probability 
Pc = P[ X P2 for the various binary systems that we considered. 

Case P3 [%] P [%] 



1 30.1 0.21 

2 29.8 0.21 

Table 5: Probabilities P3 and the total combined probability of capture and a 
circular retrograde orbit at the current epoch, P = Pi X P2 X P3, for tidal cases 
1 and 2. 

cumulative distribution of the binary's peri-Neptune distance is 
linear, the probability of disruption Pi oc rtd oc (mi + ma)"'^"', 
which is the approximate trend observed in the table. 

The value of P3 depends on the tidal models that we used. 
Its values are listed in Table|5]and the unit is again percent. Both 
cases 1 and 2 are considered. The value of P - P i x P2 x Pj, 
is also listed, where we used the average values of Pi and P2 
from the previous table. 

The total probability of a successful capture, Pc - P\ x P2, 
is on average approximately 1:140. If we then restrict ourselves 
to Triton being retrograde and circular at the current epoch, the 
probabiUty becomes P= 1:500. We can use these values to con- 
strain the number of primordial binaries, A^^, having one Triton- 
sized member. The probability of having at least one Triton 
capture if we have Nb Triton binaries is 1 minus the probabil- 
ity of its negation, that is, the probability of having no capture 
after all A^^ cases. Since the events are independent this total 
probability becomes the product of each individual probabil- 
ity of not having a capture. Each individual probability of not 
having a capture is 1 - P = 139/140 so that the probability of 
having no capture after Nb trials is (139/140)'^''. The probabil- 
ity of having at least one capture is of course its complement 
Pj„( ^ 1 _ (1 _ P)^ ^ 1 _ (139/140)''''. Since Neptune's sys- 
tem appears to be disrupted and the Uranian system appears to 
be regular, the maximum value of Ptot is approximately 50%. 
Solving for A^^ we obtain A^^ = 97 for the number of binaries 
at the planetary instability epoch with at least one Triton- sized 
member. Is this consistent with current theory? Morbidelli et 
al. (2009) claim there were Nt ~l 000 Tritons in the trans- 
Neptunian disc at the planetary instability epoch, within factor 
of a few. This implies a primordial binary population with at 
least one Triton-sized member, and with the mass ratios that 
we considered, of approximately 10%. The primordial Kuiper 
Belt population is thought to have consisted of 5%-40% bina- 
ries (Burns, 2004; Noll et al., 2008; Lin et al., 2010), so that 
our value of A^^ is consistent with this estimate, though only at 



the lowest level. However, the above fraction corresponds to 
smaller size bodies than the ones we are interested in here, and 
there is some indirect evidence that the primordial binary pop- 
ulation for heavier bodies is lower (Brown et al., 2006). Thus 
we conclude that the Agnor & Hamilton (2006) binary capture 
scenario during the planetary instability is consistent with Tri- 
ton's existence. 

However, in the above argument we did not place any re- 
striction on the final semi-major axis, eccentricity and inclina- 
tion of Triton. Placing restraints on the final semi-major axis 
and inclination will decrease the final probability, P, because a 
fourth probability, P4, comes into play. It would contain infor- 
mation about the chances of having Triton end up in a specific 
range of semi-major axis and/or inclination. This introduces the 
question of choosing a suitable bin size for the final semi-major 
axis and/or inclination that could be compatible with Triton's 
current orbit. We prefer not to do that here and instead only 
quote the probability of Neptune ending up with a circular, ret- 
rograde satellite of Triton's mass. 

Given that Neptune has a large, retrograde satellite and the 
other giant planets do not, a natural question to ask is why did 
this not happen for the other giant planets. The only reason 
we can think of is that Neptune encounters more planetesimals 
than the other giant planets (see table [TJ. Even though Uranus 
encounters a similar number, the encounters occur at a greater 
speed which decreases the probability of capture. 

4. Nereid: The fly in the ointment 

In the previous section we have presented the results of our 
numerical simulations. We concluded that the capture of Triton 
through the dissociation of a binary that had a deep encounter 
with Neptune, followed by subsequent tidal evolution that left 
Triton semi-molten (McKinnon, 1984), is enough to place Tri- 
ton on its current circular, retrograde orbit. No extra ingre- 
dients, such as the collision with a hypothetical regular satel- 
lite of Neptune (Goldreich et al., 1989), the presence of a gas 
disc (McKinnon & Leith, 1995) or a debris disc left over from 
the mutual colUsions among the members of a fictitious regular 
satellite system (Cuk & Gladman, 2005) are needed (Correia, 
2009). Does that mean that we are done? 

No, because there are several issues that we did not address. 
The first of these is what would have happened to Neptune's ir- 
regular satellites if Triton were captured during the planetary 
instability, when the other irregular satellites were captured too 
(Nesvorny et al., 2007). While Triton remained on a highly- 
eccentric orbit, with a semi-major axis comparable to that of 
the other irregular satellites, it greatly perturbs the rest of this 
population (Cuk & Gladman, 2005). Nereid in particular is dif- 
ficult to keep, and Cuk & Gladman (2005) conclude that it wiU 
be lost within 0. 1 Myr. However, the initial conditions of Cuk 
& Gladman (2005) are somewhat artificial and direct compari- 
son with our post-capture orbits is difficult. Therefore we have 
performed similar simulations, in which we place Triton on an 
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Figure 1 1 : Lifetime vs. initial semi-major axis of a fictitious swarm of irregular 
satellites of Neptune that were perturbed by Triton. The initial pericentre of 
Triton was set at 7 but the semi-major axis was varied (see titles above 
panels for values). 

eccentric orbit with semi-major axis ranging from 100 to 400 
Neptune radii, pericentre distance q ~ 1 and investigated 
how its presence affected Nereid and a swarm of other irregular 
satellites. We placed 1 000 fictitious irregular satellites around 
Neptune where the initial conditions were taken from Nesvorny 
et al. (2007). The Sun was added as an external perturber It 
turned out that Nereid is always lost through collision with Nep- 
tune within 0.1 Myr, even if it was captured with an initially 
almost circular orbit. The longest stability was found when Tri- 
ton's semi-major axis was less than 100 Rn. Using the tidal 
model, Triton needs of the order of ~100 Myr to become cir- 
cular. The eccentric phase lasts for approximately Ta = a/ a 
which for highly-eccentric orbits becomes 

'"-'■ilk) (i^) 

Thus when Triton's semi-major axis is 200 Rn and q ~ 7 R^ it 
stays eccentric for approximately 3 Myr, but if a = 100 Rn it is 
not even 2 Myr This simple approximation does not take Kozai 
mechanism into account, which increases the time by approx- 
imately one order of magnitude. Thus, with this model Triton 
stays eccentric for far too long and we lose Nereid. We have 
presented the results of some numerical simulations in Fig. [TT] 
The panels plot the time a satellite is lost vs. its initial semi- 
major axis. As can be seen, even when Triton has an orbit with 
just fl = 100 Rn, Nereid (at a - 222 Rn) is lost within 0. 1 Myr, 
even when Nereid's initial eccentricity is close to zero. Finally, 
Fig. [12] shows the original semi-major axis and eccentricity of 
the irregular satellites that are lost (bullets) and that survive 
(open circles). The big filled squares indicate Neptune's current 
irregular satellites Nereid, Halimede, Sao and Laomedeia. The 
farthest two, Psamathe and Neso, are off the scale to the right. 
As one can see, Nereid is always lost, even if it was captured 
with a low eccentricity. Only satellites captured with a low ec- 
centricity and which also have a semi-major axis a > laj sur- 
vive. 



Figure 12: Original semi-major axis vs. initial eccentricity of a fictitious swarm 
of irregular satellites of Neptune that were perturbed by Triton. Bullets show 
objects that are lost, while the open circles show objects that survive. Large 
filled squares show the positions of Neptune's currently-known irregular satel- 
lites Nereid, Halimede, Sao and Laomedeia. The initial pericentre of Triton 
was set at 7 Rn but the semi-major axis was varied (see titles above panels for 
values). 

There are five ways out of this dilemma that we can think 
of. The first is that Triton could have been captured with a very 
small semi-major axis, smaller than ~ 50 Rn- In the top panel 
of Fig. [13] we plotted the cumulative distribution of the origi- 
nal semi-major axis upon capture for all objects which reached 
an eccentricity e = 10"^ within the age of the solar system. 
The solid line is for case 1, the dashed line represents case 2. 
As one can see, only 5% of the time is Triton captured with a 
semi-major axis < 50 Rn. The median value is around 200 Rn, 
approximately where Nereid is {a = 222 Rn), which is much 
lower than the median value upon capture (Fig. [Sjl. Thus, cap- 
ture at small semi-major axis is unlikely, but not impossible. 
For reference, the bottom panel contains the cumulative semi- 
major axis distribution of objects with a > 5 Rn and e > 10"^. 
Most of these are still exhibiting Kozai oscillations and have 
undergone very little tidal evolution. 

A second possible exit strategy is to argue that Nereid was 
not captured by the mechanism of Nesvorny et al. (2007) and 
instead was a regular satellite of Neptune that was scattered out- 
wards by Triton. We performed a series of numerical simula- 
tions to test this hypothesis. We placed a large number of test 
particles on circular, equatorial orbits around Neptune with a 
maximum semi-major axis of 45 Rn, which is twice as far as 
Oberon is from Uranus. We have not witnessed any regular 
satellite being scattered to Nereid's orbit by Triton. Instead, 
they all collide with Neptune within a few tens of thousands 
of years, or they are ejected by Triton when their orbits reach 
^ < 3 Rn and their orbital angular momentum is at a minimum. 
Even if Triton could place a regular satellite on Nereid's orbit, 
we are still faced with the dilemma of keeping it there after- 
wards. In addition, we have one additional argument against 
the scattering scenario. 
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Figure 13: Top panel: Cumulative distribution of semi-major axis upon capture 
of objects that reach eccentricity e = 10"^ within the age of the solar system. 
The solid line represents case 1, the dashed line is case 2. Bottom panel: the 
same as the top panel but for orbits with final a > 5 Rn and e > 10"^ . 

The measured rotation period of Nereid is 13.6 h (Grav et 
al., 2003). This has two implications: i) either Nereid was cap- 
tured on its current orbit and thus the rotation period that we 
see is a remnant from when it formed, because it is too far away 
from Neptune to be tidally despun, or ii) it used to be on a 
13.6 h orbit around Neptune and was scattered to its current or- 
bit by Triton. A period of 13.6 h corresponds to a semi-major 
axis of approximately 3 Rn, closer in than Proteus at 4.7 R^- If 
Nereid was any further from Neptune its rotation period would 
be longer than it is now. So it is extremely unlikely that Tri- 
ton scattered a synchonous Nereid from close to Neptune to its 
current orbit without messing up the rest of the system, in par- 
ticular the satellites Proteus and Galatea. Thus a captured origin 
for Nereid is the most plausible. If this capture occured around 
the same time as that of Triton, Nereid would be lost through 
collision with Neptune. 

A third scenario has been proposed by Desch & Porter (20 1 0), 
who investigated the idea of Triton having been a satellite of a 
2 Earth-mass fictitious planet called Amphitrite. When the bi- 
nary Amphitrite-Triton suffered a close approach with Neptune, 
Triton was captured in orbit around Neptune. Desch & Porter 
(2010) speculated that Amphitrite could later collide with Uranus 
to produce its axial tilt or with Neptune itself to account for 
its excess heat radiation. Assuming that Amphitrite collided 
with Neptune and that this caused the capture of Triton, De- 
sch & Porter (2010) state the probability of capturing Triton in 
this manner is approxiately 20%-40% for encounter velocities 
at Neptune's Hill sphere of less than 3 km s '. Their typical 
pericentre distance of Triton after captre is 7 R^, for original 
Amphitrite-Triton semi-major axis shorter than 40 Rn- Unfor- 
tunately Desch & Porter (2010) do not give any information 
about the typical semi-major axis of Triton after capture, so that 
it is unclear if this capture mechanism is able to capture Triton 
at a short-enough semi-major axis in order to prevent the loss 
of Nereid. In addition, for low encounter velocities and long 



initial semi-major axis of the binary, the pericentre distance of 
Triton upon capture is large and tides may not be able to circu- 
larise it within the age of the Solar System. 

A fourth scenario involves the debris disc proposed by Cuk 
& Gladman (2005). However, as we stated in the introduction, 
it is not clear whether or not Neptune's hypothetical regular 
satellites will grind themselves down to a debris disc before 
one of them collides with Triton and shatters both Triton and 
itself. In principle the satellites collide with each other because 
Triton forces their eccentricities. This forcing is inversely pro- 
portional to the semi-major axis ratio of the regular satellites 
and Triton, and Triton's eccentricity. Furthermore, if Triton's 
semi-major axis is very long it is energetically more favourable 
for the satellites to eject Triton. Thus we envision three pos- 
sible outcomes as a function of increasing original semi-major 
axis of Triton after capture. When Triton's semi-major axis is 
short, mutual collision among the satellites, which form the de- 
bris disc, is the most likely outcome because Triton induces a 
large eccentricity in the satellites. For intermediate semi-major 
axis of Triton, the forced eccentricities of the satellites are not 
large enough for them to cross each other and Triton will hit one 
of these satellites before it is ejected. For very long semi-major 
axes of Triton, it will be ejected before a collision occurs. We 
reserve investigating this scenario for future work. 

The fifth, most plausible, scenario is that Triton was cap- 
tured and circularised before the planetary instability, some- 
thing which was already suggested by Vokrouhlicky et al. (2009) 
but for different reasons. This early capture scenario would 
solve the problem of destabilising Neptune's other irregular satel- 
lites, such as Nereid, because it was unlikely that they were 
already there. Since Neptune might not have had a circumplan- 
etary gas disc but just a cooUng, rotating envelope (Ayliff'e & 
Bate, 2009), the perturbations from Triton could have disturbed 
this system enough to prevent the formation of a regular satel- 
lite system such as that of Uranus, while decreasing its orbital 
energy at the same time. Proteus' current position would be 
indicative of the minimum distance Triton reached after its cap- 
ture prior to reaching its current orbit. In addition, at this early 
stage Neptune's obliquity may have been close to 0, increas- 
ing the likelihood of Triton ending up on its current orbit (see 
Fig. Is). Neptune's low obliquity would also have changed the 
final inclination distribution, which we have plotted in Fig. [141 
As one can see, the retrograde orbits are more evenly spread 
over a larger interval. We should add that it is possible that the 
final semi-major axis and inclination distribution of captured 
objects before the planetary instability would be diff'erent than 
that presented above, because the binary encounters with Nep- 
tune would have occurred at different velocities since Neptune 
was most likely closer to the Sun, and because it was not mi- 
grating. However, we do not think that the final results would be 
qualitatively very different from what has been presented above, 
although the probabilities discussed earlier would most likely 
change. The early capture scenario also revokes the need for 
a debris disc resulting from mutual collisions among fictitious 
regular satellites of Neptune (Cuk & Gladman, 2005) or a col- 
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Figure 14: Histograms of tlie frequency of tlie final inclination with respect to 
Neptime's equator for case 2 i.e. without Neptune's obliquity. 

lision with such a fictitious satellite (Goldreich et al, 1989) and 
only requires the minimum ingredients: a binary capture and 
tidal evolution. 



5. Summary and conclusions 

We have made an attempt at determining the origin of Tri- 
ton, Neptune's large, circular, retrograde satellite. Our work is 
based on two assumptions. The first is that Triton was cap- 
tured through a binary exchange (Agnor & Hamilton, 2006) 
and that the subsequent shrinking and circularisation of the or- 
bit occurred through tides only (Correia, 2009). We performed 
numerical simulations of the migration of Neptune in the frame- 
work of the Nice model (Tsiganis et al., 2005) and recorded the 
close encounters that the planetesimals suffered with Neptune. 
The closest of these encounters were re-enacted with binaries 
of various mass ratios to determine whether the capture of Tri- 
ton would occur, and for what parameters. The re-enaction ex- 
periments yielded a set of initial orbits for Triton just after its 
capture. These post-capture orbits were integrated using the 
tidal model of Mignard (1979, 1980) and Hut (1981), as was 
done by Correia (2009). However, we added the effects of 
Kozai mechanism induced by the Sun (Kozai, 1962) and the 
effect of Neptune's figure (e.g. Kinoshita & Nakai, 1991). The 
former will cause oscillations in Triton's eccentricity and in- 
clination. Thus Triton's orbital angular momentum, which is 
constant when only tides are taken into account, is no longer 
conserved. The perturbations from the figure of Neptune will 
overrule the Kozai effect once Triton is close enough to Nep- 
tune, so that the eccentricity and inclination oscillations cease. 
We integrated the tidal equations until Triton's orbit was circu- 
lar (e < 10"^). 

We find that a binary capture and tidal evolution are suffi- 
cient to reproduce Triton's current orbit, even though the tidal 
model would predict Triton to be closer to Neptune than its cur- 
rent position at 14.3 Rn- The probability of Neptune having a 



Triton-mass satellite is 0.7%. From this we deduce there were 
~ 100 binaries in the trans-Neptunian disc with at least one 
Triton-sized member This number is consistent with theoreti- 
cal predictions, though at the low end. 

The typical time for a final orbit at 14- to become circular 
is of the order of 200 Myr, much shorter than the age of the solar 
system or the time between the formation of the giant planets 
and the planetary instability. However, the Kozai oscillations 
caused by the perturbations from the Sun increase the time that 
Triton stays in an eccentric orbit. The perturbations from an 
eccentric Triton destabilise Nereid on a time scale of less than 
0. 1 Myr, so that its existence is in contradiction with our model. 
The most plausible exit strategies are either capture of Triton at 
small semi-major axis (a rare event), or a capture before the 
planetary instability when Nereid was not there. The early cap- 
ture of Triton would remove the need for a fast circularisation 
of its orbit, decrease the possibility of a collision with an ex- 
isting regular satellite and also increase the probability of Tri- 
ton ending up on its current orbit because Neptune's obUquity 
may have been close to zero. Thus, we suggest that Triton was 
captured shortly after Neptune's formation through a binary en- 
counter and was circularised to its current orbit through tides, 
while possibly disturbing Neptune's rotating envelope (Ayliffe 
& Bate, 2009) and preventing the formation of other regular 
satellites. 
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